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Abstract 

We present the exact supersymmetric solution of Schrodinger equation with the 
Morse, Poschl- Teller and Hulthen potentials by using the Nikiforov-Uvarov method. 
Eigenfunctions and corresponding energy eigenvalues are calculated for the first six 
excited states. Results are in good agreement with the ones obtained before. 
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1 Introduction 

Supersymmetric quantum mechanics (SUSYQM) is a framework used to determine energy 
eigenvalues and eigenfunctions of quantum mechanical problems. Various methods have 
been used in their solutions [1]. One of them is the factorization method introduced by 
Schrodinger to study the hydrogen atom problem [2]. This method was developed later 
by Infeld and Hull to categorize the analytically solvable potential problems [3]. Also it 
was used to derive in the solutions of certain non- linear equations [4]. The others are 
approximation methods known as 1/N expansion [5], S — expansion [6], supersymmetric 
WKB (SWKB) [7, 8] and variational methods [9]. About two decades ago, Gendenshtein 
introduced a concept "shape invariance" to obtain the exact energy eigenvalues of the 
Schrodinger equation [12, 13]. It provides us a simple approach to the solution when the 
potentials have shape invariance property. Recently, an alternative method known as the 
Nikiforov-Uvarov method (NU-method) has been introduced for solving the Schrodinger 
equation (SE). There have been several applications of SE with some well-known potentials 
[24, 25, 26] , Dirac and Klein-Gordon equations for a Coulomb potential by using this method 
as well [27]. 

This article is organised as follows: In Sec. (2), we give a summary of SUSYQM. In Sec. 
(3), we introduce the NU-method to solve SE. In Sec. (4), solution of Morse, Poschl- Teller 
and Hulthen potentials are obtained. Numerical results are given in the tables. Finally, we 
discuss the results in the conclusion. 



2 Summary of SUSYQM 

Supersymmetric algebra allows us to write Hamiltonians as [1] 

where The supersymmetric partner potentials V±{x) in terms of the superpotential W{x) 
are given by 

V^{x) = W^±^'-^. (2) 
The superpotential has a definition 



2 



W{x) 



h 



'2m 



dx 



(3) 



where, ^^o*^^ {x) denotes the ground state wave function that satisfies the relation 



*J°^(x) = 7Vo exp 



'2m 



W{x') dx' 



(4) 



The Hamiltonian H± can also be written in terms of the bosonic operators A and A"*" 



(5) 



where 



±^=— + W{x) 



, ,. (6) 

y/2mdx 

It is remarkable result that the energy eigenvalues of H_ and i7+ are identical except 
for the ground state. In the case of unbroken supersymmetry, the ground state energy of 
the Hamiltonian if_ is zero {e'q^ — [14]. In the factorization of the Hamiltonian, the 
Eqs. (1), (5) and (6) are used respectively. Hence, we obtain 



= {AtA-,) + Ef\ 



(7) 



Comparing each side of the Eq. (7) term by term, we get the Riccati equation for the 
superpotential Wi {x) 



Let us now construct the supersymmetric partner Hamiltonian H2 as 



(8) 



H2{x) = 



and Riccati equation takes 



r d^ , 
-2^d^^+^'^''^ 

{A^At)+Ei'\ 



(9) 



(10) 
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Similarly, one can write in general the Riccati equation and Hamiltonians by iteration as 



{AtAt)+E^\ (11) 



and 



where 



AtA-+Ef\ r»= 1,2,3,... (12) 



'2m dx dx 

Because of the SUSY unbroken case, the partner Hamiltonians satisfy the following expres- 
sions [9, 14] 

= with 4°^ = 0, n = 0,1,2,... (14) 

and also the wave function with the same eigenvalue can be written as [14] 

^(1) = A^^, (15) 

with 

= (16) 
This procedure is known as the hierarchy of Hamiltonians. 

3 The Nikiforov-Uvarov Method 

The NU-method reduces the second order differential equations (ODEs) to the hypergeo- 
metric type with an appropriate coordinate transformation x — x{s) as 
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where a{s) and a{s) are polynomials with at most second degree, and f(s) is a polynomial 
with at most first degree [24, 25, 26, 27]. If we take the following factorization 

= (18) 

the Eq. (17) becomes [27] 

a{s) y"{s) + t{s) y'{s) + A y{s) = 0, (19) 

where 

a(.)=7r(5)^(ln0(5)), (20) 



and 



Also, A is defined as 



t(s) =f(s) + 27r(s). (21) 



A„ + nr' + M^^^^ = 0, n = 0,l,2,... (22) 

The energy eigenvalues can be calculated from the above equation. We first have to 

determine 7r(s) and A by defining 

A; = A-7r'(s). (23) 
Solving the quadratic equation for 7r(s) with the Eq. (23), we get 



a' — T \ . I f a' — f^'^ 



^(5) = (^^-J ± y 1^^- j -a-rka. (24) 

Here, 7r(s) is a polynomial with the parameter s and prime factors denote the differentials 
at first degree. The determination of k is the essential point in the calculation of 7r(s). 
It is simply defined as by setting the discriminant of the square root must be zero [27]. 
Therefore, we obtain a general quadratic equation for k. 

The determination of the wave function is now in order. We consider the Eq. (20) and 
the Rodrigues relation 
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where Cn is normalizable constant and the weight function p{s) satisfy the following relation 

^^[a{s)p{s)\^T{s)p{s). (26) 

The Eq. (25) refers to the classical orthogonal polynomials that have many important 
properties especially orthogonality relation can be defined as 

rb 

/ Vnis) Vmis) p{s) = 0, Ul ^ U. (27) 
Ja 

4 Calculations 

We will use the Nikiforov-Uvarov method by getting the hypergeometric or confluent 
hypergeometric form of the Schrodinger equation with the following potentials for s — states 
only. 

4.1 Morse Potential 

The Morse potential is 

VM{r) = D « ^ - 2e-'^ ^) , (28) 

with X = {r — rQ)/rQ [23]. Here, D denotes the dissociation energy parameter and Tq is the 
equilibrium distance between nucleus. Its supersymmetric form becomes [1] 

W{x) = A - Be"" ^ (29) 
Thus, we flrst get the superpartner potentials as 

V±{x, A, B)=A^ + B^e-^ " ^ - 2B f A ^ 4^ ) (30) 



/8m, 

Prom the Eq. (1), Schrodinger equation can be written as [8] 

= (31) 

or explicitly 



6 



Using the Eq. (30) for V^, we get 



2777 r - - - n 

h a 



(33) 



where E — — , B — and C = 2S A H — -^= . By introducing a transformation 



8m , 



-s = e 



the Eq. (33) takes the form 



*"(s) + -*'(s) + 4 f/^' - - a's^l *(s) = 0, 



where 



2mB 



/32^^(£;'<0) and 7^ 



2mC 



Comparing the Eq. (35) with the Eq. (17), we obtain 

a{s) — s, f{s) — 1 and a{s) — (^—(3^ — — a^s^^ 
Substituting these polynomials into the Eq. (24), we get 



(34) 



(35) 



(36) 



(37) 



7r(s) = i^aV + (A; + 7^)5 + /32. 
The constant k is determined as 



(38) 



/ci, 2 = -7' ± 2a/9, 



and we have 



(39) 



7l{s) 



±{as-(3), for k = --f^ - 2a(3 
±{as + P), for fc = -72 + 2a/?. 
A proper value for 7r(s) is chosen, so that the function 



(40) 



r{s) = (1 + 2/3) - 2q;s, 



(41) 
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has a negative derivative [27]. By using the Eq. (22), we can find 



A = -7^ - 2aP - a 
— 2an. 



(42) 



Thus, we simply get the energy eigenvalues as 



2fe2 



E, 



n, t=0 



ah 
2m 



1 ^2 
D-in+-) 



(43) 



where D — 7^/20;. By setting a — 1, this equation reduces to the Eq. (10) as [28] for 
s — states. The Eq. (43) can be seperable. Its square term refers to the anharmonic 
oscillator correction and the other one corresponds to the harmonic oscillator solution. 

Now, we are going to determine the eigenfunctions for this potential. By considering 
the Eq. (18) and using the Eq. (20), we obtain 



By using the Eqs. (26) and (25), we obtain 



(44) 



(45) 



where p{s) — s^^ e The Eq. (45) stands for the associated Laguerre polynomials, that 
is 



where t — 2/3. Hence, we can write the wave function in the final form 



(46) 



(47) 



with s — — e" ^. It is normalizable. Using the Eq. (27), the normahzation constant can be 
found as 



Cri 



nl 



\ (n + (3 + l)(n + 2(3)V 



n 



0,1,2 



(48) 
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4.2 Poschl-Teller Potential 



The Poschl-Teller potential is 



Vpt{x) = 



Uo 



cosh^ax^ 

where Uq — \ {\ — 1) > Q Also, its superpotential potential is [1] 



W{x) = Ata.nh.ax. 
Prom Eq. (2), we get its superpartners 



V±{x,A) = A^-A{At 



Thus, we can write the Schrodinger equation as 



ah 



'2m ) cosh'^ax 



Till \ 2m 



where E ^E^- and A \ A^ 



E- 



ah ^ 
J2m , 



K 



cosh^ax. 



^!{x) = 0, 



Introducing a transformation 



s — tanh ax, 



we rewrite the Eq. (52) 



-2s 



S2) 



where 



By comparing the Eq. (54) with the Eq. (17), we determine polynomials as 



a{s) = {l-s^), f{s) = -2s and a{s) = -p^ + j'{l - s^). 
Substituting them into the Eq. (24), we obtain 



7r{s) = ±i//52-72(l-s2) + A;(l-s2). 



The constant k is determined in the same way. Therefore, we get 



±/3, for k — j^ 

7r(s) = < 

^ ±(3s, for A; = 72 - (3^. 
Here we choose the proper value, so that 

has a negative derivative. Prom the Eq. (22), we calculate 



(58) 



(59) 



A = 7^-/?2-/? 
= n'^ + n + 2n(3. 

Hence, the energy eigenvalues are found as 



Er,. ^A'- 



h'^a^ [ , 1. 1 



2m 



-(n+-) + -0^ 



(60) 



The wave function ^'„(a;) is obtained from the Eq. (18) by taking tt{s) 
follows. We first get 



(61) 
-/3s as 



0(s) = (1 - sY/', 
and using the Eqs. (26) and (25), we find 



(62) 



Vnis) 



p{s) ds^^ 



;i - srp{s) 



(63) 



where p{s) = (1 — s^Y . The Eq. (63) stands for the Jacobi polynomials as 



Vn 



(64) 



Hence, ^n(2^) can be written in the following form 



^nix)=Cr.il-Sy/' Pi^^^\s), (65) 

with s — tanha x. Considering the Eq. (27), the normalization constant is obtained as 



1 /(2n + 2/3+1) 



2/3 (n + /3)! 



n!(2n + /3)!. 



(66) 



where n, (3 > 0. 
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4.3 Hulthen Potential 

This potential can be solved exactly for s — states only. This is due to the similarity between 
Coulomb and Hulthen potentials. It plays an important role in the applications of quantum 
scattering theory. The Hulthen potential is given by [23] 

< - -^0(^4^, (67) 

where S = 1/a, is the screening parameter. We get the supersymmetric form the potential 
for s — states [10] 

^'-~^ + \i-l-s.y (68) 

Here, a and b are arbitrary constants. We can also write the supersymmetric partner of the 
potential as [11] 

The second term in Eq. (69) behaves hke centrifugal barrier [23]. The SE has the form 



2m 

Using the transformation 



ATH r 1 

*"(x) + ^\E- Vf{x)] ^{x) = 0. (70) 



s = e-' ^ (71) 

we rewrite 



where 



-{e' + /3')s' + {2e' + /3' - Y)s - ^(s) = 0, (72) 



2 2mE — 2 2mVo 2 2mVQ 



By comparing the Eq. (72) with the Eq. (17), we get 



(7(s) =s(l-s), f(s) = l-s and a{s) ^ -{e^ + P')s^ + {2e^ + - Y)s - e\ (TA) 
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Substituting them into the Eq. (24), we obtain 



n{s) = -^s ± ^^4 (^£2 + /32 - A; + i ) s2 - 4 {2e^ + - -f^ - k) s + Ae^. (75) 
From Eq. (75) k is determined as 



k,,2^-l' + P^±ey/l + 4^. 
Following the same procedure, we get 



(76) 



1 1 



(2e + ^1 + 472^ s-2s , 
and the energy eigenvalues for the supersymmetric Hulthen potential becomes 



(77) 



E. 



n, e=o 



-Vn 



2n P 



n= 1,2, 



Here, /32 = 2Vo/5^ with (fi = m = 1) and n = 



(78) 



. If the limit 7 — > is 



chosen, the energy eigenvalue reduces to the form obtained from the usual solution of the 
Hulthen potential. 

The wave functions can now be obtained similarly from the Eq. (18). Using the Eq. 
(20), we have 



where = 1 + VI + 472. The Eqs. (26) and (25) lead to 

Here, p{s) = s^^ (1 — sY~^. It stands for the Jacobi polynomials as [25] 



(79) 



(80) 



y„(.)c^p(^^' ^-^)(l-2s). 



(81) 



Thus, the final form of the wave function can also be written in terms of the Jacobi poly- 
nomials resulting 



^n{.x) = Cn (1 - 5)^/2 p{2e, /.-l)^^ _ 3^)^ 



(82) 
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with s — e , and also the normahzation constant C„. 

5 Conclusions 

We have obtained the exact supersymmetric solution of some central confining potentials 
by applying the Nikiforov-Uvarov Method. The eigenfunctions and corresponding energy 
egenvalues of the these three well-known shape invariant potentials, i.e. Morse, Poschl- 
Teller and Hulthen are calculated analytically. All the wave functions are physical. We 
present numerical results in tabular form for i = 0. In Table I, we list energy eigenvalues of 
H2 molecule by taking D = 4.7446 eV, a = 1.9425 A'^ and m = 0.50391 amu for ground 
state and first eleven excited states. The energy difference between the successive states 
decreases. In Table II, the six excited energy states of the the potential are given for various 
values of n and A with A = {1 + y^l + 4 A (A - l))/2. In Table III, the first five excited 
energy levels are tabulated for s — states with different values of screening parameter S. 
Our results are in good agreement with the ones obtained by the other methods. 
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Table I: Vibrational (non-rotating) energy eigenvalues of the Morse potential for H2 molecule. 



Quantum no (n) 


Eigenvalues [28] 


En (Our work) 


En+l — E„ 





-4.47600 


-4.47610 


0.51354 


1 




-3.96256 


0.48226 


2 




-3.48030 


0.45100 


3 




-3.02930 


0.41968 


4 




-2.60962 


0.38842 


5 


-2.22050 


-2.22120 


0.35713 


6 


-1.86330 


-1.86407 


0.32587 


7 


-1.53740 


-1 53820 


29456 


8 




-1.24364 


0.26330 


9 




-0.98034 


0.23201 


10 




-0.74833 


0.20073 


11 




-0.54760 
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Table II: Eigenvalues of the Poschl- Teller potential with^ = 
setting by = a = 2m = 1, ^2 = and k = 7^) in Eq. (61) [29]. 



(1 + VI + 4 A (A 



\ 


— El — E2 


— E3 


-E4 


-E5 


—Eg 


2.0 


0.0 1.0 


4.0 


9.0 


16.0 


25.0 


6.0 


0.0 


1.0 


4.0 


9.0 


16.0 


12.0 




0.0 


1.0 


4.0 


9.0 


20.0 






0.0 


1.0 


4.0 


30.0 








0.0 


1.0 


42.0 










0.0 
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Table III: Eigenvalues of the Hulthen potential for several values of screening parameter S, 
substituting Vq = 6/ A in of the Eq. (78). 



Oiiantiim no ( yi\ 


-E„ [301 


^exact 


— E„ ( Our workl 




6 = 


0.002 




1 


0.4990005 


A /lAAAAAr 

0.4990005 


A /lAAAAAT 

0.4990005 


o 
Z 


(J.iz4U(JzU 


A T O /I AAOA 


A 1 O /I AAOA 

U.lz4UUzU 


6 


0.0o45dU1 


0.0545601 


A A c: /I c:/:? A1 

0.0545601 


A 

4 


O.OoOzooO 


0.0302580 


0.0302580 


o 






n nni 9'inn 

U.UUIZOUU 




o — 


u.ux 




1 


0.4950125 


0.4950125 


0.4950125 


2 


0.1200500 


0.1200500 


0.1200500 


3 


0.0506681 


0.0506681 


0.0506681 


4 


0.0264501 


0.0264500 


0.0264500 








U.UluOlZO 






0.05 




1 


0.4753125 


A A <~TV ^1 CiV 

0.4753125 


0.4753125 


z 


A 1 A1 o c: AO 

0.1012503 


A 1 A1 O CT AA 

0.1012500 


A 1 A1 O C: AA 

0.1012500 


6 


0.0333746 


0.0333681 


0.0333681 


A 

4 


0.0113035 


A AIIOr'AA 

0.0112500 


A A1 1 AP'AA 

0.0112500 


cc 
o 






n nn9si 9^^ 

U.UUZOlZO 




o — 






1 






0.4900500 


2 






0.1152000 


3 


0.0460057 


0.0460056 


0.0460056 


4 


0.0220512 


0.0220500 


0.0220500 





U.UiiZ004 


n ni 1 9c:nn 

U.UiiZoUU 


U.UiiZooU 




s = 


0.2 




1 


0.4049962 


0.4050000 


0.4050000 


2 


0.0450856 


0.0450000 


0.0450000 


3 






0.0005556 


4 






0.0112500 


5 






0.0450000 
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